Phase transitions of a tethered surface model with a deficit angle term 
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Nambu-Goto model is investigated by using the canonical Monte Carlo simulations on fixed con- 
nectivity surfaces of spherical topology. Three distinct phases are found: crumpled, tubular, and 
smooth. The crumpled and the tubular phases are smoothly connected, and the tubular and the 
smooth phases are connected by a discontinuous transition. The surface in the tubular phase forms 
an oblong and one-dimensional object similar to a one-dimensional linear subspace in the Euclidean 
three-dimensional space R^. This indicates that the rotational symmetry inherent in the model 
is spontaneously broken in the tubular phase, and it is restored in the smooth and the crumpled 
phases. 

PACS numbers: 64.60.-i, 68.60.-p, 87.16.Dg 
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I. INTRODUCTION 

A considerable number of studies have been conducted 
on the phase structure of the elastic membrane model 
of Polyakov-Kleinert and Helfrich Q, S H S H S S 
H H lis [H m m El. Xhc Hamiltoman of the 
model includes not only the Gaussian term but also 
a bending energy term, which can make the surface 
smooth. Thus it has been recognized that there are a 
smooth phase and a crumpled phase in the model of 
spherical topology. Numerical studies have also been 
made to understand the phase transitions in the teth- 
ered model and in the fluid model on triangulated sur- 
faces jj III HI Hi, .19, 20, 21, 23, 24, .25, 26, 27], 
where the tethered and the fluid models are those de- 
fined on fixed connectivity surfaces and on dynamically 
triangulated surfaces, respectively. Except those smooth 
surfaces and crumpled surfaces, some linear structure 
mm US 113 and branched polymer surfaces [Hill HI 
have also been recognized in the context of surface mod- 
els. 

On the other hand, little is known about a model of 
tubular surfaces. A tubular surface is considered as one 
of the basic forms of real physical membranes [33 . US] . 

Reviewing references on tubular surfaces, we must 
recall anisotropic surface models that have been con- 
structed to understand tubular surfaces [3, • A tubu- 
lar phase is realized in the anisotropic model due to an 
anisotropic bending modulus. 

However, there has been no study that tried to under- 
stand tubular surfaces from an isotropic surface model. 
Therefore it will be interestng to study an isotropic teth- 
ered surface model corresponding to the Nambu-Goto 
string . The purpose of this study is to understand the 



phase structure of the tethered surface model of Nambu 
and Goto with a deficit angle term, which is obtained 
from the co-ordination dependent term. No bending en- 
ergy term is included in the Hamiltonian, whereas the 
deficit angle term serves as a curvature energy smoothing 
the surface. We will show numerically that the model has 
three distinct phases; crumpled, tubular, and smooth. 
The crumpled and the tubular phases are connected by 
a higher-order transition, whereas the tubular and the 
smooth phases are connected by a first-order transition. 

We consider that oblong tubular surfaces can be seen 
in the Nambu-Goto model. A discretized Hamiltonian 
of the Nambu-Goto string is given by the sum of the 
area of triangles, which corresponds to the Gaussian term 
of the Polyakov-Kleinert model. The area term of the 
Nambu-Goto model imposes a constraint only on the area 
of triangles, and hence all the triangles become oblong 
and form spiky configurations. In fact, it is well known 
that the partition function of the Nambu-Goto model is 
not well defined 0|. However, it is possible that the 
partition function of the Nambu-Goto model changes to 
a well-defined one if some additional term is included in 
the partition function, as was suggested already in Ref. 
[38| . Then, it is expected that such oblong triangles may 
form tubular surfaces in such a well-defined model if the 
additional term tends to modify the spherical surface to 
a tubular one. 



II. MODEL 



The area energy Si is defined by 
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where A/^ is the area of the triangle A in a triangulated 
surface of spherical topology. This energy 5*1 of Eq. ^ is 
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a straightforward discretization of the Nambu-Goto ac- 
tion denoted by 5* = / (PxyTg^ where g is the determinant 
of the first fundamental form on the world surface swept 
out by strings. 

The partition function Z of the Nambu-Goto surface 
model is defined by 

N 

Z{a) = / []rfX,cxp(-5), 
S{X,T)^S,-aS3, 53 = ^log((5./27r), (2) 

i 

where the Hamiltonian S is given by a linear combination 
of the Gaussian term Si and the deficit angle term 6*3, 
in which Si is the sum of the angles of vertices meeting 
at the vertex i. S{X,T) denotes that S depends on the 
embedding X and the triangulation T fixed on a uniform 
lattice, whose construction will be described further in 
Sec. mil The center of the surface is fixed to remove the 
translational zero mode. Z{a) of Eq. ||2Jl denotes that 
the model is dependent on the parameter a which is the 
coefficient of the deficit angle term 5*3. It should also be 
noted that the Hamiltonian S{X, T) is defined only with 
intrinsic variables of the surface and hence independent 
of the extrinsic geometries. 

The deficit angle term S3 has a deep connection with 
the integration measure dXi js^, |4fl Ell ■ We have to re- 
mind ourselves that dXi can be replaced by a weighted 
measure dXiqf , where qi is the co-ordination number of 
the vertex i, and a is considered to be a = 3/2. Con- 
sidering that qi is a volume weight of the vertex i, we 
think it is possible to assume a as an arbitrary number. 
Moreover, the co-ordination number qi can be replaced 
by the vertex angle 5i according to 

TlidXiq^ IljdXj exp(Q! ^ log 6i). 

i 

The constant term log 27r is included to normalize ^3 
in Eq. so that S'3 = when 5i = 2tt at every vertex. 
Thus we have the expression S3 in Eq. (01. Note that 
53 = not only on the flat surface but also on cylindrical 
surfaces. 

It should also be noted that the deficit angle term S3 
is expected to play a nontrivial role in the discrete model 
which has a finite number of vertices N. We expect that 
53 is not irrelevant in surface models such as the Nambu- 
Goto model or the Polyakov-Kleinert model. In fact, it 
was reported by Ref. '2511 where the phase transition of 
the fluid model of Polyakov and Kleinert depends on the 
coordination dependent term '^i^ogqi. 

The unit of physical quantities can be explained as 
follows: The Hamiltonian S in Eq. ^ is obtained from 
S = aSi—aS3 by assuming the surface tension coefflcient 
a as a = 1. It should be noted that the choice a — 1 
represents not only a redefinition of a as a/ a but also a 
choice of the unit of length as i/fcT/a = 1 , where T is 
the temperature, and k the Boltzmann constant. Thus 



the unit of a becomes kT/a. The choice of y^kT/a = 1 
for the unit of the length is possible because of the scale 
invariant property of the partition function in Eq. (jSJ . 

III. MONTE CARLO TECHNIQUE 

The triangulated surfaces, on which the model is de- 
fined, are uniform in the co-ordination number q. On 
the uniform lattice, the number of vertices Nq of g = 5 
is 7V5 = 12, and all other vertices are of q = 6. These 
lattices were obtained by Monte Carlo (MC) simula- 
tions with the dynamical triangulation for a model whose 
Hamiltonian is defined by Sg + ^5*2 — ctS3 (g) with suf- 
ficiently large a, where Sq = J2i ^ is the Gaussian 
term, >S'2 = X)i(l ~ cos^i) the bending energy term, and 
'^'3(9) = X]iioS9i' which is different from S3 in Eq. |(5J. 
The well defincdness of the uniform lattice must be con- 
firmed, since there are finitely many uniform lattices con- 
structed in this technique for each N . We have first con- 
firmed that a fixed connectivity surface model, which is 
defined by SG + bS2, is well defined on such uniform lat- 
tices. In fact, the specific heat for the bending energy 52 
is independent of the choice of the uniform lattice. 

The canonical Metropolis technique is used to up- 
date X. The position Xi is moved to a new posi- 
tion Xl = Xi + AXi , where AXi is randomly chosen 
in a small sphere. X'^ is accepted with the probability 
Min[l,exp(-Z\5)], where AS = S'(new) - S'(old). The 
radius of the small sphere for AXi is chosen at the be- 
ginning of the simulations to maintain 35 ~ 60 % of ac- 
ceptance rate; almost all MC simulations are done on 
about 50 % of acceptance rate. 

A lower bound 10^^ Aq is imposed on the area of tri- 
angles in the update of X, where Aq is the mean area 
of the triangles computed at every 250 MCSs (Monte 
Carlo sweeps) and ^0 is constant due to the relation 
Si/N = 1.5. However, the areas are almost free from 
such constraint, because the areas of almost all triangles 
are larger than 10~^^o throughout the MC simulations. 
No constraint is imposed on the bond length. 

IV. RESULTS 

We first show in Figs, da) and (b) the mean square 
size X^ defined by 

i i 

where X is the center of the surface. X^ in Fig. QJa) 
represents that the size of surfaces continuously increases 
with increasing a, and also represents that the shape of 
the surfaces rapidly changes at a ~ 100. On the other 
hand, shown in Fig. ^h) clearly represents some dis- 
continuous transition, where X^ abruptly changes. The 
dashed lines drawn vertically on the data of the iV= 1000 
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FIG. 1; (a) X"^ obtained at relatively small a in the vicinity of 
the boundary between the crumpled and the tubular phases, 
(b) X'^ obtained at relatively large a in the vicinity of the 
boundary between the tubular and the smooth phases. The 
unit of X^ and a is kT/a, where a is the surface tension. 



surface in both of the figures represent the phase bound- 
aries, on which we focused our attention in this paper. 

The convergence speed of MC is very low in the tubu- 
lar phase close to the smooth phase. The total number 
of MCS at a = 20000 on the iV = 1500 surface is about 
2 X 10^, where 1.5 x 10^ MCSs were discarded for the 
thermalization. This is the reason why we use surfaces 
of size up to = 1500. The reason of the low conver- 
gence speed seems due to a straight-line structure of the 
surface, which will be shown below. Since the vertices 
can move only along the line, the surface deforms very 
slowly. On the contrary, the convergence both in the 
smooth phase and in the crumpled phase is relatively 
faster than that in the tubular phase. In the simulation 
the expected relation Si/N = 3/2 is satisfied in the con- 
figurations reached after the thermalization at every a. 




500 1000 N 500 1000 n 



FIG. 2: (a) Log-log plot of vs iV at a = 50 (crumpled 
phase), and = 200 (tubular phase), and (b) those obtained at 
a close to the phase boundary of the discontinuous transition. 
The unit of X^ is kT/a. 

Figure ^a) shows log-log plots of X"^ against N ob- 
tained at a = 50 and a = 200. Plots of X^ against N 
in Fig. |2tb) denoted by tubular {smooth) were obtained 
below (above) the discontinuous transition point in each 
iV as shown previously in Fig. H^b). The straight lines 



plotted in Figs. Efa) and (b) are those fitted by 

oc N^/", (4) 

where H is the Hausdorff dimension. From the slope of 
the plotted lines, we have 

iJso = 7.24 ± 0.48, H^mo = 1-93 ± 0.01, (5) 

where H^q and i?smo were obtained from the data de- 
noted by a = 50 in Fig. Ufa) and those by smooth in Fig. 
|2Ib), respectively. Those results are in agreement with 
our expectation. In fact, H is expected to be very large in 
the crumpled phase, and it is also expected to be 7f = 2 in 
the smooth phase. Moreover, we have 772oo = 1-80 ± 0.02 
and Htuh = 1.22±0.03, which were obtained from the data 
denoted by a = 200 in Fig. [2Ia) and those by tubular in 
Fig. m^b), respectively. H200 and -fftub slightly deviate 
from H = 2 which is confirmed in the case of branched 
polymer surfaces |26| , where surfaces randomly stretches 
and hence are rotationally symmetric. 




FIG. 3: Snapshots of A'^ = 1000 surfaces obtained at (a) a = 50 
(crumpled), (b) a = 200 (tubular), (c) a = 12000 (tubular), 
and (d) « = 16000 (smooth). Surfaces in (a), (b), and (d) are 
drawn in the same scale, which is different from that in (c). 

Snapshots of = 1000 surfaces are shown in Figs. 
ini;a)-(d) obtained at a = 50, a = 200, a = 12000, and 
a = 16000. Figures Ola), (b), and (d) were drawn in 
the same scale, which is different from that in Fig. ^c). 
The axis direction of the surface in Fig. |2tc) is spon- 
taneously chosen. The direction of the axis remains al- 
most unchanged throughout the MC simulation. Thus 
we find no tubular surface bending in the tubular phase 
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for N < 1500. The straight-line structure shown in Fig. 
OIc) is expected to survive even at a sufficiently large N. 
The reason is because both 5*3 and oblong triangles tend 
to straighten the surface. 
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FIG. 5: Maximum bond length L 
a, and at (b) relatively large a. The unit of L 
and that of a is kT/a 
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max is sJkTja, 



FIG. 4: Normalized distribution hiJJ) of the bond length L 
obtained at (a) a = 50 (crumpled), and a — 200 (tubular), 
and at (b) a = 12000 (tubular), and a = 16000 (smooth), on 
iV = 1000 surfaces. The unit of L is y'kT/a. 

Figures 0{a) and (b) are normalized distribution h{L) 
of the bond length L sampled at every 500 MCSs in the fi- 
nal 2x10'' MCSs on A^= 1000 surfaces. The normalization 
of h{L) is given by '£i h{L)/AL = 1, where AL = 0.02, 
and the sum runs over all bonds and hence ^ - 1 be- 
comes identical with Nb the total number of bonds. The 
dashed and solid curves denoted by crumpled and tubu- 
lar in Fig. Efa) were obtained at a = 50 and a — 200 
respectively, and those denoted by smooth and tubular in 
Fig. mb) were obtained at a = 16000 and a = 12000, 
respectively. 

We note that h{L) obtained on surfaces of size other 
than 1000 are exactly identical with h{L) in Fig. 0{a) 
if a is identical with each other. Moreover, h{L) obtained 
at the smooth phase, denoted by smooth in Fig. 0Jb), is 
independent of both N and a in the smooth phase. While 
h[L) is dependent on a in the tubular phase close to the 
crumpled phase, it is almost independent of both N and 
a in the tubular phase close to the smooth phase. The 
fluctuations of the surface size in the tubular phase close 
to the smooth phase are relatively larger than those in 
the tubular phase close to the crumpled phase. 

It should be noted that the distribution h{A) of the 
area A of the triangles, which is not presented in a figure, 
is universal in a sense that h{A) is independent not only 
of N but also of a. h{A) is not influenced even by the 
discontinuous transition. In fact, h{A) obtained in the 
tubular phase is exactly identical not only with that in 
the smooth phase but also with that in the crumpled 
phase. 

It must be checked that the size of triangles is negligi- 
ble compared to the size of surfaces at sufficiently large 
N . In order to see that the maximum bond length L^ax 
is considerably smaller than the size of surfaces in the 
tubular phase, we plot in Figs. IHJa) and (b) L,„ax ob- 
tained in the final 2 x 10^ MCSs on each surface. We 
find in Fig. EJa) that L,„ax continuously increases with a 



at the boundary between the crumpled and the tubular 
phases, and that Lmax at each a is almost independent 
of N . In the tubular phase at a = 200, Lmax is smaller 
than the surface length L^: 

Ls -41 [N = 600, a = 200), 
Ls -56 [N ^ 1000, a = 200). (6) 

The length Ls were obtained by Ls — \J {Lfj, where (L^) 
was obtained in the tubular phase close to the crumpled 
phase. While Lmax in Fig. EJa) is almost independent 
of N , Ls increases with iV at a = 200 as shown in Eq. 
©. Hence it is expected that ^max/is ^ in the limit 
^ oo at least in the tubular phase close to the crumpled 
phase. 

FigureE^b) shows that Lmax in the tubular phase grad- 
ually increases as N increases. However, we find that 
^max is considerably smaller than Lg in the tubular phase 
close to the smooth phase. In fact, we have Lmax — 43 for 
AT = 1000 and Lmax — 54 for = 1500, which are smaller 
than the length of the surfaces shown below: 

Ls ~ 122 {N = 1000, a = 12000), 
Ls~ 172 (iV = 1500, a = 20000). 

Therefore it is also expected that Lmax/^s ^ in the 
limit N—^oo everywhere in the tubular phase. 

Finally, we plot in Figs. [HJa) and (b) the bending 
energy 5*2 = ^ cos 9i), which reflects a smooth- 

ness of surfaces and it is not included in the Hamilto- 
nian. While S2/NB continuously changes against a in 
Fig. Ela), it is clearly discontinuous in Fig. EJb). These 
results indicate that the tubular phase is smoothly con- 
nected to the crumpled phase and discontinuously con- 
nected to the smooth phase. The higher-order nature 
of the transition between the tubular and the crumpled 
phases has also been seen in the specific heat Cgg = 
(a^/N) {(Si)-(S3f). In fact, although Cs^, has a peak 
at a ~ 100, there was no growth of the peak with increas- 
ing A^. 

The bending energy S2/NB is not an order parameter 
regarding to the tubular phase because it is also nonzero 
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FIG. 6: The bending energy 5'2/iVs = Ei(l " cos ei)/NB at 
(a) relatively small a, and at (b) relatively large a. Nb is the 
total number of bond. The unit of a is kT/a. 



in the crumpled phase. However, S2/NB plays a role of 
order parameter, as can be seen in Fig. Etb), as far as we 
confine ourselves to the transition between the tubular 
phase and the smooth phase. 



phases, and a higher-order transition between the tubular 
and the crumpled phases. The surface forms an oblong 
and one-dimensional object in the tubular phase. It is re- 
markable that the rotational symmetry or the symmetry 
of isotropy inherent in the model is spontaneously broken 
in the tubular phase. 

An important point to emphasize is that both terms 
area Si and deficit angle ^3 are the cause of such variety 
of phases. Moreover, it is quite likely that the straight- 
line structure in the tubular phase survives even at suffi- 
ciently large N, because not only the oblong triangles but 
also the deficit angle term can make the surface tubular. 

Further numerical studies on the fluid model and on 
the model with extrinsic curvature would give us hints 
to clarify the phase diagram of the Nambu-Goto surface 
model with the deficit angle term. 
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V. SUMMARY AND CONCLUSION 

We have investigated the phase structure of a teth- 
ered surface model of Nambu-Goto embedded in R^, and 
found that there are three distinct phases: smooth, tubu- 
lar, and crumpled. Moreover, the model undergoes a 
first-order transition between the smooth and the tubular 
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